Flux state and anomalous quantum Hall effect in square Kondo lattice 
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The anomalous Hall effect (AHE) around the flux state in square Kondo lattice is investigated. By 
introducing the lattice distortion and local chirality, the square Kondo lattice can break the parity 
symmetry and time reversal symmetry spontaneously, and thus generate a topological nontriviality 
in the band structure associated with the AHE. Moreover, a possible realization of this AHE in 
multiferroic TbMnOs is discussed. 

PACS numbers: 75.10. Jm,75.30.Mb,75.47.Lx 



In condensed matters, novel spin orders often lead to 
novel physical phenomena. For instance, some multifer- 
roics have a spiral spin order which breaks the spatial 
inversion symmetry and gives rise to ferroelectric polar- 
ization, which's origin is completely different from con- 
ventional ferroelectricity [H, 0] ■ This spiral spin order can 
be scaled by a vector spin chirality (VSC) Si x Sj and 
has become an important concept in the physics of spin 
current and spin liquid 0, Besides the VSC, there is 
another scalar spin chirality defined by: 
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which breaks the parity (P) and time reversal (T) sym- 
metry and was first proposed by Wen et al 5] . It is clear 
that chirality x would be nonzero for a noncoplanar spin 
order. 

On the other hand, recent studies showed that the non- 
coplanar spin order is relevant to the intrinsic anomalous 
Hall effect (AHE) observed in Kondo lattice system, in 
which the noncoplanar background spin texture acts as 
a gauge field for itinerant electrons propagating in the 
lattice [i, 0, H, S, 113, [m, E, O, Q, \M- Generally 
speaking, the intrinsic AHE has the topological origin 
and can be characterized by the Berry phase and Chern 
number. To manifest this mechanism, the system should 
break the P-symmetry and T-symmetry spontaneously 
and simultaneously 3, lB| • For a Kondo lattice model, 
the T-symmetry can be violated for those spin configura- 
tions with the local spin chirality, while the P-symmetry 
can be broken in some geometrically frustrated lattices. 
For instance, Ohgushi et al once discussed the AHE on 
Kagome lattice where a finite local spin chirality in the 
three-site unit cell can generate nonzero Hall conductance 
0. This model has been extended to other geometri- 
cally frustrated systems such as the three-dimensional 
(3D) pyrochlore lattice 0, B 0|- Therefore, a geomet- 
rically frustrated lattice with spin chirality would be of 
significance in terms of AHE physics. 

Unfortunately, the square lattice usually has no geo- 
metrically frustrated structure, and thus it is not easy 
to violate the P-symmetry since the chiralities on ad- 
jacent plaquettes tend to cancel each other due to the 



lattice symmetry and thus the AHE becomes hard to re- 
alize in the square lattice. However, the square lattice 
(and its distorted forms) takes up the majority in the 
practical Kondo lattice materials such as colossal mag- 
netroresistance manganites and Fe-based pnictide su- 
perconductors Besides, the realization of AHE on 
the square lattice is also fundamentally important and 
physically interesting 



In this paper, the AHE on the Kondo square lattice 
is realized theoretically by introducing two mechanisms 
to break the P-symmetry. One is to induce some lattice 
distortions which can lead to the change of hopping am- 
plitude of itinerant electrons. The other is to construct 
a special unit cell which breaks the P-symmetry. These 
two mechanisms are different from the previous consid- 
ered spin-orbit interaction (SOI) which can directly gen- 
erate topological nontrivality in the band structure and 
associate with the AHE [l^, [11, [2l| . Our discussion will 
be primarily restricted to a 'flux' state at the half fill- 
ing of one-band Kondo model. We hope that these two 
mechanisms can be alternating approaches (other than 
the SOI mechanism) to the AHE in the square lattice, 
and eventually realized in some real materials. 

The Hamiltonian of one-band Kondo lattice model on 
the two-dimensional (2D) square lattice can be written 
as: 113 



H 



NN NNN 



-Ji ^ Si • Sj 



NN 



J2 

NNN 



(2) 



where the first term describes the electron hopping be- 
tween the nearest-neighbor (NN) sites, and the hopping 
amplitude ti is taken as the energy unit. The second 
term is the electron hopping between the next-nearest- 
neighbor (NNN) sites. In the following, t2 is arbitrarily 
set as 0.25 as an example, since the AHE result is qual- 
itatively independent of its exact value as long as it is 
nonzero, which will be further dicussed below. The third 
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term is the Hund coupling linking the itinerant electrons 
with the background spins S (assumed classical and nor- 
malized as |S| — 1) where Jh is the coupling factor. The 
last two terms are the antifcrromagnetic (AFM) superex- 
changes between the background spins with J\ and J2 as 
the coefficients for the NN and NNN sites respectively. 
This model has been extensively investigated for various 
transitional metals' oxides, and more details of this model 
can be found in Ref . [17 . 

For the third term, by applying a canonical transfor- 
mation, the Hamiltonian can be simplified by adopting 
the site-dependent spin-polarization axis. In the Jh — > 00 
limit, the spin of the hopping electron is forced to align 
parallel to the on-site S, and thus the hopping terms 
in the Hamiltonian can be transferred into the form 
t'^j'^clcj, with the effective hopping integral: 
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where t can be ti or ^2- ^ and ip are the polar coordinates 
of spin S. The phase factor ipij can be viewed as the 
gauge vector potential and 0ij is the angle between 



and S," 



When the itinerant electrons move along a 



closed loop, they can feel the induced gauge flux which 
is indistinguishable from the magnetic flux. This gauge 
flux is related to the spin chirality and leads to the AHE 

For this model, the groud state can be calculated 
with the variational method, namely, by comparing the 
ground state energies of several preset spin configura- 
tions . Around half filling, when the Ji=J2 > 0.15, the 
ground state is the fiux state. This state has four sites in 
the unit cell, as shown in Fig. 1(a). In one plaquette, the 
neighboring background spins are perpendicular to each 
other. When an itinerant electron travels around the pla- 
quette, it can acquire an additional tt flux. In fact, this 
flux state was reported earlier in some similar systems 
(2^ [23! . The fiux phase can have many degenerate states 
by shifting and rotating the spin structures. Thus, a spe- 
cific state, with (Si, S3) along the z-axis and (S2, S4) 
along the a;-axis, will be adopted in the following study, 
which forms a spin order in the x-z plane. In practical 
calculation, to remove the degeneration and stabilize this 
state, a small anisotropic energy H2 = ^aS^^ (sub- 
script a=l-4, Ki=K:i=-Qm, K2=Ki=Qm) is also con- 
sidered, which favors (Si, S3) along the z-axis and (S2, 
S4) in the x-y plane. In fact, the fiux state is with copla- 
nar spin order and the AHE conductance is forbidden. 
Therefore, additional contributions should be included 
for allowing the AHE. 

In real materials, due to the ionic size mismatch 
or competing exchange interactions, the square lattice 
would be distorted more or less. For instance, in multi- 
ferroic RMn03, the Mn-O-Mn chain is distorted in non- 




FIG. 1: (Color online) (a) The flux state with the four-site 
unit cell, (b) The lattice is distorted with atoms 1 and 3 
moving in reversed direction with the new positions shown as 
open circles. The NN hopping amplitude t\ is varied corre- 
spondingly. The arrows on bonds indicate the signs of the 
phases of the t^j^ ■ (c) The band structure of Eq. (5) with the 
parameters sets tri = 1, tr2 = 0.25, d = 1 and <j> ~ 7r/6. (d) 
Gauge flux density of the 3-th band of (c). 

centrosymmetric manners caused by the Dzyaloshinsky- 
Moriya (DM) interaction, giving rise to the staggered 
Mn-O-Mn angles P, [2^ . For simplicity, a simple lattice 
distortion mode is adopted, as shown in Fig. 1(b), where 
the first and third cations are displaced from the origi- 
nal positions along the opposite directions. Due to the 
x-y symmetry, the cation displacement along x and y are 
assumed to be the same. In a first order approximation, 
the corresponding hopping amplitude varies linearly with 
the tiny lattice distortion, and thus the NN hopping am- 
plitude t[ = ti±d where d is the tiny amendment caused 
by the distortion. The infiuence to the NNN hopping 
amplitude t2 is not considered because the exact value 
of NNN hopping is not qualitatively important to obtain 
the AHE . Take one chain of the lattice for example, the 
hopping amplitude becomes staggeringly ordered and the 
P-symmetry of the lattice is broken. 

With this P-symmetry broken configuration, the spin 
order may be not exactly confined on the x-z plane, i.e., 
spins Si and S3 tilt slightly from the z-axis and are not 
parallel to each other. Therefore, the local spin order 
may become non-coplanar, with each plaquette having 
a gauge fiux penetrating it. In this condition, the NN 
hopping can be approximately modified to the following 
form [6|: 



t' 



eff 



trie 



(4) 



where t^i is the rcnormalizcd NN hopping amplitude. 
We take tpij — (j){^<j)) for the hopping direction along 
(opposite to) the arrow direction as shown in Fig. 2(b). 
In addition, the renormalized NNN hopping amplitude is 
set as a real constant tr2 for simplicity. The Hamiltonian 
matrix for this model can be written in the momentum 
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space: 
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H(k) = 
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(5) 



where fi{k)=tri cos{kx/2)+idsm{kx/2), 

f2{k)=triCOs{ky/2)+idsmiky/2), /3(fc)=tr2[cos((fca; - 
ky)/2)+cos{{kx + ky)/2)]. Now, the AHE conductance 
can be calculated. The contribution to the AHE 
conductance from each band is written as 0, [l^ : 



— / d'^kz ■ Vk X A„(fc) 
Jbz 



h 2TTi 
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where A„(k) =< nk\V\nk > is the vector potential de- 
fined with the n-th wave function, B„(k) is the gauge 
fiux density and C is the so-called first Chern number. 
At d = case, namely, the ideal square lattice without 
any distortion, the P-symmetry is maintained, leading to 
zero Chern number for each band. At d 7^ case, the 
P-symmetry is broken. The calculation indicates that 
the Chern number for each band is C=[0, 0, 1, — 1] at 
(f) < 7r/4, and C=[l, -1,0,0] at </> > 7r/4 (at = 7r/4, 
the T-symmetry is conserved, corresponding to the fiux 
state). For instance, at ^ = 7r/6, the band structure is 
shown in Fig. 1(c), and the gauge flux density B(k) of 
the 3-th band is shown in Fig. 1(d). 

For Hamiltonian Eq. (2), the nonzero NNN hopping 
term is essentially important to generate the topological 
nontriviality in the band structure. This can be intu- 
itively understood as follows. For Eq. (5), in the large 
d-limit, fi{k) ~ idsin(kx/2X under the unitary trans- 
formation UHU^ H' [2J|, the Hamiltonian matrix 
Eq. (5) can be further decoupled into the form: 



' h2{k) 



where /ii(fc) and h2{k) are both 2x2 matrixes. 
For the matrix hi{k), we have: 
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where /4(fc) = sin(fca;/2) -|-sin(fcj^/2), f^ik) — sin(/c^/2) — 
sin(fcy/2). Around k = (7r,7r), the electron can be con- 
sidered as a generalized Dirac fermion and the effective 
Hamiltonian hi{k) {h2{k) can be treated in a similar way) 
is 

h^{k) = -^fc^fc' + ^^(fc'2 _ k'^)av + 2dsin0fT^ 

(9) 



where k'^ = kx — tt, k'y = ky — it. The general form 
of Eq. (10) was thoroughly addressed in Ref. [H, and 
the corresponding Chern number for the upper and lower 
bands is C = i:2sgn{tr2/d) (at (j) = 7r/6) where sgn is 
the sign operator. If ^^2 = 0, C = for both bands, 
indicating that the NNN hopping term is indispensable 
in generating the Hall conductance. 

Consequently, one can argue that the lattice distortion 
provides an effective method to break the P-symmetry, 
and yet the non-coplanar spin order as the ground state 
is absolutely necessary. In fact, the non-coplanar spin 
order can also be accomplished by further taking into ac- 
count the frustrated magnetic interaction. For instance, 
by adding the third neighbor (3rd N) superexchange in- 
teraction i?3 = X^srdAT "^3^« ' ^3 distorted Kondo 
lattice model described by Eq. (2), the original flux state 
at the half filling evolves into a state in Fig. 2(a). In 
this condition, the unit cell expands to eight-site and is 
formed with two interlaced square sublattices, compared 
with the original four-site unit cell (formed by Se, S3, S7 
and S4). For this spin configuration, we can straightly 
calculate the Hall conductivity at zero temperature by 
using the Kubo formula [l^ : 



e 

h 2ni 



-E 



< nk\Jx\mk >< mk\Jy\nk > —h.c. 



[£n(fc) 



(10) 

where \nk > is the occupied state and \mk > is the empty 
state, e denotes the eigen-energy. Jx (Jy) is the x (y) 
components of the current operator. The calculation is 
performed in momentum space, with the summation over 
all the eight bands. For example, with a finite J3 (see Fig. 
2(a)'s caption), (T2;j,=0.13 e^/h at half filling, clearly in- 
dicating that this ground state does exhibit the AHE. 
Besides, the Hall conductance as the function of con- 
duction electron density n is shown in Fig. 2(b), with 
the spin configuration fixed as Fig. 2(a). The curve in 
Fig. 2(b) fluctuates dramatically, showing that the AHE 
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is sensitive to n 

Interestingly, additional investigation of the ground 
state in Fig. 2(a) indicates that besides the lattice dis- 
tortion, this special unit cell structure also breaks the 
P-symmetry spontaneously. Thus for this unit cell struc- 
ture, with the local spin chirality, the band structure 
can lead to nonzero Chern number and Hall conductance 
even without any lattice distortion or NNN hopping. For 
the local spin order shown in Fig. 2(c), the eight bands 
are topologically nontrivial and the corresponding Chern 
number C = [1,-1,0,0,0,0,-1,1], with the gauge flux 
density of the 1-th band shown in Fig. 2(d). In fact, due 
to the frustrated magnetic interaction, the unit cell could 
involve more sites and become even larger with the parity 
violation satisfied. For instance, for Eq. (2), with large 
Ji {J2=Ji) and low n, the variational calculation indi- 
cates the ground state is with the coplanar spiral order. 
This coplanar spiral order can evolve into the 3-D spi- 
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distortion and the other is the locally nontrivial spin or- 
der. Both of these mechanisms break the P-symmetry 
and can generate the AHE spontaneously. 

The authors acknowledge with E. Dagotto, Q. H. 
Wang, N. Nagaosa and R. Shindou for fruitful discus- 
sions. 



FIG. 2: (Color online) (a) The ground state with the parame- 
ter sets Ji = J2 = 0.04, J3 = 0.025, ti = 1, f2 = 0.25, d = 0.4 
and n = 0.5. The unit cell expands to eight-site here, (b) The 
Hall conductance as a function of the conduction electron den- 
sity n with the spin order shown in (a), (c) An example of 
the nontrivial spin order in the special eight-site unit cell, (d) 
Gauge flux density B(/c) of the I-th band with the spin order 
shown (c). 

ral order by including the H2 and H3 term. However, in 
this case, the unit cell also expands dramatically, and the 
band calculation becomes much more challenging. The 
ground state we address here seems to be a special case 
with a relatively small unit cell. 

Although the calculation of Eq. (2) and the addressed 
two mechanisms are more or less theoretically oriented, 
their realization in real systems is still possible. For ex- 
ample, multiferroic TbMnOa is a promising candidate to 
illustrate these two mechanisms and observe the AHE. 
In TbMnOa, besides the GeFeOs-type distortion, which 
forms the Mn-O-Mn zigzag chain, the DM interaction is 
also present 25|. In this case, due to the spiral order of 
the background spin, the Mn-O-Mn zigzag chain is dis- 
torted with all the oxygen ions displaced in the same di- 
rection [25j. Therefore the adjacent Mn-O-Mn angles are 
different, leading to the staggered order of the effective 
Mn-Mn hopping amplitude and broken parity symmetry. 
For TbMnOs, neutron scattering experiments confirmed 
that the background i2g electrons form a coplanar spiral 
order. To excite the noncoplanar spin order, a magnetic 
field with its direction different from the spin order plane 
can by applied, allowing the nonzero Hall conductance. 
In addition, the ferroelectric polarization in TbMnOa 
aligns along the direction of Gy x (S^ x Sj). Therefore, 
the noncoUinear ferroelectric polarization should be ex- 
pected for this non-coplanar spin order. Moreover, the 
Tb3+ cations can be doped by other -1-4 cations, which 
can modulate the itinerant electrons density and directly 
control the AHE f2l'|. However, for these calculations, a 
more practical two-orbital model should be employed. 

In conclusion, we have studied the Kondo lattice model 
on a square lattice with frustrated super-exchange inter- 
actions. The calculated nonzero Hall conductance is at- 
tributed to two distinct mechanisms, one is the lattice 
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